clear all
close all

data = spiral_dataset(0:0.01:15);
errors = [abs(data(:, 3)) abs(data(:, 1))+abs(data(:, 2)) ones(size(data, 1), 1)]*1;
for i = 1:size(data, 1)
    data(i, :) = gsamp(data(i, :), diag(errors(i, :)), 1);
end

plot3(data(:, 1), data(:, 2), data(:, 3), '.')
axis equal
